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Abstract. The serendipitous discovery of an unidentified extended TeV y-ray source close to the galactic plane named 
HESS J1303-631 at a significance of 21 standard deviations is reported. The observations were performed between February 
and June 2004 with the H.E.S.S. stereoscopic system of Cherenkov telescopes in Namibia. HESS J1303-631 was discovered 
roughly 0.6° north of the binary system PSRB1259-63/SS 2883, the target object of the initial observation campaign which 
was also detected at TeV energies in the same field of view. HESS J1303-631 is extended with a width of an assumed intrinsic 
Gaussian emission profile of cr = (0. 16 + 0.02) ° and the integral flux above 380 GeV is compatible with constant emission over 
the entire observational period of (17 ± 3) % of the Crab Nebula flux. The measured energy spectrum can be described by a 
power-law dN/dE ~ E~ Y with a photon index of T = 2.44 ± 0.05 sta t ± 0.2 syst . Up to now, no counterpart at other wavelengths is 
identified. Various possible TeV production scenarios are discussed. 



16 

17 
18 



Key words. Gamma rays: observations - Galaxy: disk 



1. Introduction 



3. Dataset 



Cherenkov telescopes of ground based y-ray astronomy 
achieve a very high sensitivity. Due to the limited field of view 
the pointing of the telescopes has in general to be decided 
on the basis of observations in other wavelengths. However, 
the comparatively large field of view of instruments such as 
HEGRA (4.3° full angle) and H.E.S.S. (5° full angle) together 
with the stereoscopic observation mode allows nevertheless a 
scan of a part of the sky (Aharonian et al. 2002b, 2005b) and 
in particular allows a search for unknown sources in the field 
of view of individual pointings. In this way the first unidenti- 
fied TeV y-ray source TEV J2032+4130 has been discovered 
by HEGRA (Aharonian et al. I2002al I2004at in archival data 
and now, with the H.E.S.S. telescopes, the second unidenti- 
fied source HESS J 1303-631 was discovered, followed by more 
unidentified sources found in the H.E.S.S. galactic plane scan 
(Aharonian et al. 2005b}. The H.E.S.S. telescopes have origi- 
nally been directed to search for TeV y-ray emission from the 
binary system PSRB1259-63/SS2883 near periastron begin- 
ning in February 2004. This binary system has been detected 
at TeV energies (Beilicke et al. [2004) which is the subject of 
a parallel paper (see Aharonian et al. 2004c). Surprisingly, an- 
other TeV y-ray source located at a position roughly 0.6° north 
of the position of the binary system was discovered in the same 
field of view, given the name of HESS J 1303-631 (see Fig. FT1. 
The detection and basic features of this new source as well 
as a search for possible counterparts in other wavelengths and 
discussion of possible TeV y-ray production scenarios are re- 
ported in this paper. 



2. The H.E.S.S. Cherenkov telescopes 

The H.E.S.S. (High Energy Stereoscopic System) collabora- 
tion operates an array of four imaging atmospheric Cherenkov 
telescopes (IACTs) optimised for an energy range of y-rays be- 
tween 100 GeV and 20 TeV. The telescopes are located in the 
Khomas Highland in Namibia (23°16'18"S, 16°30'1"E) at a 
height of 1 800 m above sea level. Each telescope has a 107 m 2 
tessellated mirror surface (Bernlohr et al. 2003; Cornils et al. 
2003) and is equipped with a 960 photomultiplier tube cam- 
era with a field of view diameter of ~ 5° (Vincent et al. 2003) 
which allows searches and studies of TeV y-ray sources in sky 
regions of more than 3° x 3° per pointing. The telescopes are 
operated in a coincident mode (Funk et al. 2004 1 assuring that 
an event is always recorded by at least two of the four tele- 
scopes allowing for stereoscopic reconstruction of the shower 
parameters. More information about H.E.S.S. can be found in 
Hinton (20041. 



Send offprint requests to: Matthias Beilicke, 
e-mail: matthias . beilicke@desy . de 

* Universite Libre de Bruxelles, Faculte des Sciences, Campus de 
la Plaine, CP230, Boulevard du Triomphe, 1050 Bruxelles, Belgium 
** now at Department of Physics and Astronomy, Univ. of Sheffield, 
The Hicks Building, Hounsfield Road, Sheffield S3 7RH, U.K. 
*** UMR 7164 (CNRS, Universite Paris VII, CEA, Observatoire de 
Paris) 



The data were taken between February and June 2004 with the 
fully operational H.E.S.S. IACT array. The average zenith an- 
gle of the observations was 42.7° yielding an energy thresh- 
old of Ethr = 380 GeV, defined by the peak y-ray detection 
rate of a Crab-like spectrum after event selection cuts. The 
observations were performed in the wobble mode, tracking a 
position shifted by +0.5° in Declination or Right Ascension 
with respect to the nominal source position (in this case the 
PSR B 1259-63/SS 2883 position of the initial observation cam- 
paign), allowing for an unbiased simultaneous background de- 
termination for not too extended sources. Following the detec- 
tion of HESS J1303-631 in the PSR B1259-63/SS 2883 field of 
view, the telescope tracking positions were changed to new 
sky positions in May 2004, optimised for both sources by 
choosing the wobble modes with respect to a position located 
between PSR B1259-63/SS 2883 and HESS J1303-631. In the 
right hand panel of Fig.^the initial pointing positions (filled 
stars) as well as the new pointing positions (empty triangles) 
are shown. 

The data were selected by standard quality criteria (stable 
weather and detector status) leaving 54.5 hours of data (48.6 h 
detector life time) for the final analysis. For the data taken be- 
tween February 26 and March 5, 2004 (7.6 hours), one of the 
four telescopes was excluded from the analysis due to technical 
reasons with the event read-out of the telescope. The remaining 
data (March to June 2004) were analyzed using the full array of 
four telescopes. The raw data were subject to the standard cal- 
ibration (Aharonian et al. 2004b) and Hillas parameter-based 
analysis (see Aharonian et al. 2005a for more details). 



4. Analysis and Results 

4.1. Discovery of TeV y-ray emission from 
HESS J1 303-631 

Shortly after the discovery of TeV y-ray emission from the bi- 
nary system PSRB1259-63/SS2883 the highly significant ex- 
cess of the second source HESS J1303-631 was found in the 
data. The discovery sky map resolving both TeV sources is 
shown in the left hand panel of Fig. ^ The number of ex- 
cess events was calculated using the ring background model 
in which the background is determined from a ring region with 
a radius r > 0.5° centred around the putative source position. 
Both TeV y-ray source regions were excluded from the back- 
ground region to avoid background contamination by real y-ray 
events. In order to improve the angular resolution for the inves- 
tigation of the position and extension, for the sky map at least 
three images instead of two were required per event. 

4.2. Consistency checks 

Since with the discovery of HESS J 1303-631 for the first time 
in the H.E.S.S. data and in TeV y-ray astronomy in general, 
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Fig.l. Left: The sky map showing both TeV y-ray sources: HESS J1303-631 and PSRB1259-63/SS 2883. The Galactic plane 
is also indicated. The number of events are integrated within a circle of < 0.14° for each of the correlated bins. Right: The 
function which was fitted to the (uncorrected) excess sky map is drawn as contour lines. Also shown are the different H.E.S.S. 
tracking positions: Filled stars indicate the initial tracking positions (February - May) and the empty triangles the optimised 
positions which were used since May, 2004. 



two sources were found in the same field of view 1 a wide range 
of consistency checks was applied and successfully passed un- 
derlining the celestial origin of the measured excess. For this 
purpose the data were divided into different subsets and the 
results (source strength, excess position, etc.) were compared. 
The subsets were chosen according to the different observation 
periods (from February to June, 2004) and the different tele- 
scope tracking positions (compare Fig[2 right). Different cuts 
on the telescope multiplicity m (m > 2, m = 2, m > 3, m = 3 
and m = 4) for individual events as well as different cuts on 
the image amplitude were applied to the data. Also, individual 
telescopes in turn were excluded from the analysis. In all these 
checks the excess behaved as a genuine y-ray source; more de- 
tails can be found in Beilicke et al. (I2005i . 

Finally, an important check concerns the distribution of the 
mean reduced scaled width (MRSW) parameter which is used 
for the y-hadron separation (see Aharonian et al. 2005a 1. The 
MRSW is obtained by averaging the reduced scaled widths of 
the individual images of an event. The reduced scaled width 
is defined as the difference (measured in standard deviations) 
between the measured and the expected width of an image. 
The expected width is obtained from Monte Carlo simula- 
tions for y-ray showers, being parameterized as a function of 
the image amplitude, the shower impact parameter, and the 
zenith angle. The distribution of the MRSW parameter of the 
HESS J1303-631 excess is shown in Fig. |2 together with the 
distributions for y-ray events from simulations and for hadronic 
events from a background measurement. The MRSW distribu- 
tion of HESS J 1303-631 is compatible with y-ray events. 
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Fig. 2. Distribution of the mean reduced scaled width (MRSW) 
parameter (see text) for the HESS J1303-631 excess (data 
points). For comparison, the distributions for simulated y-ray 
events (dashed line) and for hadronic events (histogram) are 
also shown. The excess measured from HESS J1303-631 is 
compatible with y-ray events. 



1 Meanwhile, two sources in one field of view have also been ob- 
served in other pointings with H.E.S.S., e.g. in the galactic center re- 
gion and the galactic plane scan. 



All results obtained from the different investigations are 
in good agreement. They were cross-checked with the tem- 
plate background model (Rowell 2003 1 and were also con- 
firmed by an alternative model analysis technique (de Naurois 
et al. 2003 1. The celestial origin of the measured excess from 
HESS J 1303-631 can therefore be considered as established. 
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Fig. 3. Distribution of ON-source events (solid histogram) 
and normalized OFF-source events (filled histogram) vs. 
the squared angular distance A0 2 between the reconstructed 
shower direction and the object position. The expected ON- 
distribution for a point-like source (TeV point spread function, 
PSF) is indicated by the dashed line; it is normalized to give 
the same number of excess events as HESS J 1303-631. 



4.3. Excess and sky position 

To obtain the position of HESS J1303-631 a function describ- 
ing the excess of the two sources in the field of view was fitted 
to the uncorrected sky map which was generated with a cut on 
at least three images per event in order to obtain a better angular 
resolution. The excess of HESS J 1303-631 was fitted by a 2D 
elliptical function (two-dimensional Gaussian function with a 
cr for the Right Ascension and one for the Declination as well 
as a covariance). Simultaneously, the PSRB1259-63/SS2883 
position was fitted by the H.E.S.S. point spread function (PSF) 
for TeV y-rays to avoid a systematic influence on the fit by ex- 
cess events from PSRB1259-63 leaking into the HESS J1303- 
63 1 region. The PSF can be described by a 2D double Gaussian 
function a ■ exp(^p-) + b ■ exp(^i) which was derived from 
Monte Carlo simulations. The fit range covers a region of 
2° x 2° centred at the HESS J1303-631 position. The param- 
eters of the PSF (crj, cr 2 and a lb) which was applied for 
the position of PSRB1259/SS2883 were left free. The con- 
tour lines of the function are shown in the right hand panel 
of Fig. [U The position of HESS J1303-631 was found to be 
a = 13 h 03 m O:4 + 4 S .4 and 6 = -63° 11 '55" ±31" (J2000.0). 
The length and width of the fitted ellipse are within errors the 
same and clearly larger than that expected for a point source 
(see below). Therefore, the HESS J1303-631 excess is com- 
patible with an extended and rotationally symmetric struc- 
ture. The position of PSRB1259-63/SS2883 obtained from 
the same fit was found to be a l259 = 13 h 02 m 49 s .2 + 3*6 and 
<5i259 = -63°50'274 + 21'.' 1 which is in good agreement with 
its nominal position of ai259, n om = 13 h 02 m 47 s .7 and £i259,nom = 
-63°50'8'.'8. The width of the PSRB1259-63/SS2883 excess 
was found to be compatible with the width of a point-source 
(Aharonian et al. I2004ct . Fixing the position of PSRB1259- 
63/SS2883 in the fit to its nominal sky coordinates does not 
change the results for the HESS J 1303-631 position and ex- 
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Fig. 4. The differential energy spectrum of HESS J 1303-631. 
A fit to the data with a power law (solid line) is shown. The 
power-law corresponding to a spectrum of 1 % of the Crab 
Nebula is also indicated (dashed line). 

tension within statistics. The systematic pointing uncertainty 
of the H.E.S.S. telescopes is estimated to be ~ 20" for Right 
Ascension and Declination. 

The distribution of the number of events in squared angular 
distance A0 2 measured between the reconstructed shower di- 
rection and the derived HESS J1303-631 position is shown in 
Fig-El The PSF of the H.E.S.S. detector for TeV y-rays is also 
shown (dotted line). A single Gaussian function describing the 
intrinsic source profile of HESS J 1303-631 was folded with the 
PSF. The folded function was fitted to the excess distribution 
resulting in an intrinsic width of <t hessji3 ( )3 = (0.16 + 0.02) 
with a^ 2 /d.o.f. = 21/42. 

Adjusting the angular cut to the derived extension (A0 2 « 

A0 pomts rc + <4essji303 ~ 0.05 deg 2 ) one obtains for the 
HESS J 1303-631 position 2469+1 19 excess events correspond- 
ing to a significance of 21 standard deviations calculated fol- 
lowing Li & Ma fE%3l 

4.4. Energy spectrum and light curve 

An energy spectrum was derived using the angular cut of 
A© 2 = 0.05 deg 2 . To obtain the full flux integrated over the 
whole emission region of HESS J 1303-631 the flux normalisa- 
tion was corrected for the derived source extension (assuming 
an intrinsic Gaussian emission profile). Since this correction 
depends on the exact - and possibly energy dependent - shape 
of the source, it introduces a systematic error. Together with 
the effects of atmospheric extinction variations and energy cal- 
ibration of the detector the systematic error on the flux was 
estimated to be in the order of ~ 30 %. 

The spectrum is shown in Fig.0] It was fitted by a power- 
law dN/dE = No ■ (.E/lTeV)~ r with a resulting photon in- 
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Fig. 5. The HESS J 1303-631 light curve covering February un- 
til June, 2004. Shown is the integral flux <£(£ > 380 GeV) vs. 
time in night-by-night bins. The fit of a constant function (solid 
line) results in a^ 2 /d.o.f. of 35/35 being compatible with con- 
stant emission. 



0.2 syst and a normalization of 



dex of T = 2.44 + 0.05 stat 

N = (4.3 + 0.3 sta t) • lO-^cm^s-" 1 TeV" 1 with a ^ 2 /d.o.f. of 
27/9. The integral flux above 380 GeV was calculated to be 
0(£ > 380 GeV) = (1.2±0.2 stat )-10- n ctrT 2 s -1 corresponding 
to (17 + 3) % of the flux of the Crab Nebula. Spectra obtained 
from the different tracking positions as well as from the dif- 
ferent observation periods were found to be compatible within 
statistical errors. A fit of a curved spectral shape, e.g. with a 
power-law and a cut-off, yields a better x 2 - However, at the 
present understanding of the systematic effects for extended 
sources observed with asymmetric wobble positions (comp. 
Fig-E right), we prefer not to quote numbers at this stage but to 
leave it to more detailed studies on morphology and spatially- 
resolved energy spectra which are underway. 

A light curve of the integral flux <$>(E > 380 GeV) was de- 
rived on a night-by-night basis and is shown in Fig.|5] To obtain 
the integral flux, the count rates for each night were corrected 
using the effective areas and an assumed power-law as obtained 
from the overall differential energy spectrum. A fit of a con- 
stant function to the light curve results in a^ 2 /d.o.f. of 35/35 
and therefore indicates constant emission from HESS J1303- 
63 1 during February until June, 2004. 



5. Search for possible counterparts 

5.1. Search in catalogues 

As HESS J1303-631 is a new TeV source, a search for pos- 
sible counterparts at other wavelengths was performed, using 
the VizieR service (Ochsenbein et al. 2000 1 at the Centre de 
Donnees astronomiques de Strasbourg. The closest EGRET 
(> 100 MeV) source, 3EG J 1308-61 12, has its centroid at an 
angular distance of 2.1° (RA = 197°18, Dec = -61°22), with a 
nominal 95 % confidence radius of 0.71° (Hartman et al. 1999 1. 
This source is flagged as possibly extended or multiple and 
confused, but examination of the actual position confidence 
contours confirms that it can be ruled out as a counterpart. 
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Fig. 6. 843 MHz SUMSS radio map of the region around 
HESSJ1303-631 (marked as the large circle of radius 0.16°). 
Also marked are the positions of the X-ray sources (WGA and 
RXP1-5) as well as the radio sources (PMN1-5, Gl and G2 for 
Hn regions and PSR1-5 for pulsars) listed in Tabled Circles or 
ellipses surrounding individual sources indicate the catalogued 
source extension (convolved with the relevant radio beam size). 

Table 1. X-ray (ROSAT) and radio sources at a distance of less 
than 30' from the centre of HESS J 1 303-63 1 . The radio sources 
marked PMN are from the Parkes-MIT-NRAO catalogue, those 
marked G are H n regions, and those marked PSR are from the 
ATNF pulsar catalogue. Shown are the source name, the ID 
used in Fig. |6j the distance to the HESS J1303-631 centroid 
and the J2000.0 source coordinates. 



source name 


ID 


dist. 


RA 


Dec 


X-ray sources (ROSAT) 








1WGAJ1303. 8-6319 


WGA 


9.1' 


13 03 50.9 


-631906 


2RXPJ130257.9-632750 


RXP1 


15.9' 


13 0257.9 


-6327 50 


2RXPJ130540.5-631839 


RXP2 


19.2' 


13 05 40.5 


-63 18 39 


2RXPJ130622.7-632210 


RXP3 


24.9' 


13 06 22.7 


-63 2210 


2RXP J 130200.4-633730 


RXP4 


26.4' 


13 02 00.4 


-63 37 30 


2RXPJ130337. 8-634045 


RXP5 


29.1' 


1303 37.8 


-63 4045 


radio sources 










PMN J 1302-6257 


PMN1 


14.3' 


1302 38.9 


-625748 


PMN J 1300-6253 


PMN2 


24.6' 


13 00 34.4 


-625343 


PMN J 1303-6339 


PMN3 


28.2' 


130352.6 


-63 39 29 


PMN J 1259-6333 


PMN4 


29.7' 


1259 59.1 


-63 33 37 


PMN J 1302-6321 


PMN5 


9.4' 


1302 56.7 


-632120 


H ii regions 










G 303.9-0.1 


Gl 


25.4' 


125956.6 


-62 57 22 


G 303.9-0.7 


G2 


28.9' 


130007.3 


-63 3320 


pulsars 










PSR J 1302-63 13 


PSR1 


4.9' 


130219.2 


-631329 


PSR J 1303-6305 


PSR2 


6.9' 


13 03 00.0 


-63 05 01 


PSRJ1301-6310 


PSR3 


10.4' 


130128.3 


-631041 


PSRJ1301-6305 


PSR4 


10.6' 


130145.8 


-63 05 34 


PSRJ1305-6256 


PSR5 


22.6' 


13 05 28.0 


-63 56 39 



In the X-ray domain, the only catalogued sources found 
with VizieR within 30' of the centroid of HESS J1303-631 are 
from the Second ROSAT (0.24 - 2.0 keV) Source Catalogue of 
Pointed PSPC Observations (5 sources) or its WGACAT ver- 
sion (1 additional source). These ROSAT sources are listed in 
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Tabled and their positions are marked in Fig. [5] The concentra- 
tion of these sources to the South of HESS J 1303-631 reflects 
the non-uniform sensitivity of PSPC observations across this 
field. Dedicated X-ray observations will be necessary to make 
a more complete search for counterparts. 

As HESS J 1303-631 is within the Galactic plane, optical or 
infrared sources are too numerous and have not been examined. 

To search for counterparts in the radio domain, the 
843 MHz map of the region around HESS J1303-631 was ob- 
tained from SUMSS survey (Bock et al.[W99 1 and is shown as 
a grey-scale image in Fig. |6] Those relatively bright sources 
within 30' of HESS J1303-631 which are also identified in 
the 4.85 GHz PMN (Parkes-MIT-NRAO) survey catalogue 
(Wright et al. 19941 are marked as circles or ellipses 2 in Fig. [6] 
and are listed in Table □ along with PMN J 1302-6321 which 
is the only catalogue source within the 1-cr source extension 
of HESS J 1303-631, but is not detected in the 843 MHz map. 
The region overlapping the centre of HESS J 1303-631 appears 
devoid of conspicuous radio sources. 

The nearest catalogued supernova remnant, Kes 17 
(G304. 6+0.1) , lies at a distance of 36' from HESS J1303-631 
(Green 2004 1. VizieR also finds two Hn regions within 
a 30' radius, G 303.9-0.1 and G 303.9-0.7 (Paladini et al. 
2003 1. Finally, version 1.21 of the ATNF (Australia Telescope 
National Facility) Pulsar Catalogue 3 (Manchester et al. 2005 1 
lists five pulsars within a 30' radius around HESS J 1303-631. 
The most compelling positional coincidences of radio sources 
with HESS J1303-631 appear to be with some of the pulsars, 
discussed below. 

5.2. Possible pulsar wind nebula associations 

Pulsars can have associated very high energy (VHE) y-ray 
emission from the nebulae of high-energy particles (mostly 
electrons and positrons) which they produce, as the well-known 
example of the Crab Nebula shows. The TeV radiation is 
thought to arise from inverse Compton (IC) scattering of lower- 
energy ambient photons by these energetic electrons, which 
also emit synchrotron radiation at lower energies, typically be- 
tween the infrared and X-ray bands. In the case of a relatively 
young (950 year old) pulsar like the Crab, the pulsar wind neb- 
ula (PWN) is fairly compact and appears point-like in VHE 
y-rays, but in older pulsars the emission might be considerably 
more extended, due to the expansion of the PWN itself (e.g. van 
der Swaluw et al. 2001 1 or due to the escape of the high-energy 
particles out of the synchrotron-emitting PWN, as proposed by 
Aharonian et al. (11997b . 

To examine the possibility that HESS J 1303-631 might be 
associated with such an extended nebula from one of the posi- 
tionally compatible pulsars, the latters' energetics are now con- 
sidered. The total available energy output of a pulsar is given by 
its spindown luminosity E. It is believed that the major fraction 
of this energy is ultimately converted to accelerated particles, 

2 The objects marked PMN1 and PMN2 are consistent with point 
sources, and their diameter reflects the PMN beam size, while PMN3 
and PMN4 are flagged as extended. 

3 http : //www . atnf . csiro . au/research/pulsar/psrcat 



Table 2. Parameters of the pulsars possibly associated with 
HESS J 1303-631. Listed are the pulsar name, its period P, es- 
timated distance D, spindown age t and flux F sc / (see text). 
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Fig. 7. Hi image (~ 35° x 20°) at v = -24 km s _I of the shell 
GSH 305+01-24 (McClure-Griffiths et al. EfflHV The member 
stars of Cen OBI, location of HESS J1303-631 and approxi- 
mate outer boundaries of GSH 305+01-24 are indicated. 

and the energetic electrons and positrons can radiate their en- 
ergy away efficiently in the form of synchrotron and inverse 
Compton emission. In such a scenario the measured VHE en- 
ergy flux from a source at the distance D can be compared to 
^sd = E/(4nD 2 ), whereby the conversion efficiency is expected 
to be much less than 100%, i.e. the measured VHE flux is ex- 
pected to be much less than F s( j. 

Table |2] lists the relevant parameters for the five pulsars 
included in Table [2 as derived from version 1.21 of the 
ATNF Pulsar Catalogue. For comparison, the energy flux of 
HESS J1303-631 between 0.3 and lOTeV implied by its best- 
fit spectrum is 2.1 x 10~ n ergcirT 2 s -1 . Some caveats should 
be kept in mind when comparing this number with the listed 
values of F^: uncertainties in the estimation of D for individ- 
ual pulsars can be significant; the comparison also implicitly 
assumes that the y-ray emission is isotropic, but this is almost 
certainly a valid assumption for an extended source such as 
HESS J 1303-631. Finally, delayed energy release of particles 
injected in the past, when the spindown luminosity was higher, 
could in principle increase the apparent efficiency; but this is 
unlikely to be a large effect in most evolutionary scenarios. 

With this in mind, it can be seen that most of the pul- 
sars listed in Table |3 can be ruled out as possible counter- 
parts for HESS J1303-631 due to insufficient energetics. Only 
PSR J 130 1-6305 has enough energy to power the whole of the 
VHE emission, and even in this case the conversion efficiency 
would have to be unusually high. 
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5.3. Molecular clouds, energetic stellar winds from 
Cen OB1, and SNR association 

If HESS J 1303-631 is a result of cosmic ray interactions with 
ambient gas, one can draw the following conclusions: the fact 
that the spectrum of this source is significantly harder com- 
pared to the E~ 2J diffuse (local) galactic cosmic ray spectrum, 
indicates that one does not see an interaction of this galactic 
cosmic ray component with regions of enhanced ambient tar- 
get gas density, but rather of a source of cosmic rays, which 
is located within or near HESS J 1303-631. If such a source is 
also in contact with a dense molecular cloud, the increased tar- 
get density results in an increased rate for the production of 
VHE y-rays as discussed by Aharonian et al. ( 1990 1. 

Molecular clouds can provide the conditions for both Type 
I and Type II supernova events to occur, with the former re- 
sulting from accelerated accretion induced collapse of white 
dwarfs due to the increased gas density, and the latter result- 
ing from the evolution of massive stars in star formation re- 
gions towards supernovae. Another possibility is that stellar 
winds from OB stellar associations, interacting with molecular 
clouds, may also be sources of y-ray production, as suggested 
for TeV J2032+4130 by Aharonian et al. J2002b> . or, such stel- 
lar winds may provide the required injection for further accel- 
eration in a nearby SNR shell from the same or another nearby 
OB association (Montmerle 1979 1, called supernova-OB asso- 
ciations (or SNOBs). There is a possible correlation between 
SNOBs and unidentified y-ray sources as found by Montmerle 
( fT979l in the COS-B data, and in the EGRET data by Kaaret, 
& Cottam ( fT996l . 

As a first step the nature of molecular clouds in the direction 
of HESS J 1303-631 is investigated, as reflected by CO transi- 
tion line measurements, which gives (i) a kinematic distance 
based on galactic rotation curves (sometimes one will get more 
than one solution), and, (ii) the column density of gas corre- 
sponding to the region from where the CO emissions are seen. 

It is possible to fit the observed spectral shape of 
HESS J 1303-631 with a hadronic model. A hadron spectral 
index of 2.0 combined with an exponential cut-off at 20 to 
50 TeV would even account for the slight curvature apparent in 
the measured spectrum. To reproduce the observed TeV flux, 
one has to assume a total cosmic ray energy of Ecr — 2 to 
5 x 10 49 d^ v J n\ erg, where d\^ c = d/l kpc is the distance and 
n\ = n/l cirT 3 is the associated target density. 

The procedure of deriving n from the velocity integrated 
CO temperatures (W(CO)) = J Tdv is well-known, provided 
that a scale size for the cloud is specified (see e.g. Dame et al. 
(1993) and references therein). The deep CO survey of molec- 
ular clouds in the southern Milky Way by Bronfman et al. 
(1989) covers the 300° to 348° longitude and -2° to +2° galac- 
tic latitude ranges, with a spatial resolution of 0.125°, which is 
smaller than the resolved size of HESS J 1303-631. 

Extracting CO temperatures T centred at / = 304.24° and 
b = -0.36° (the centroid for HESS J1303-631) from this sur- 
vey, results in three distinct spatial coincidences in velocity 
space. The merits of these will be discussed below. 

(i) It is unlikely that the Coalsack (v = -4 to -2kms _1 ), 
at a distance of ~ 175 pc (Nyman et al.[l989), is a candidate 



environment for HESS J1303-631 based on two aspects: there 
is no evidence of star formation in this dark cloud (Nyman et 
al.|T989 1, and a bright counterpart should have been detectable 
at longer wavelengths, given the proximity of this cloud. 

(ii) For the second CO solution the reader is refered to 
Figure (11) of Bronfman ( 1989 >, who shows a clump in the 
velocity interval ~ -30 to -15kms _1 along the line-of-sight 
cutting through the centroid of the TeV source. The tempera- 
ture was integrated along this velocity interval to give W(CO), 
showing a plateau at the location of HESS J 1303-631. This 
quantity, if averaged over the surface of the resolved TeV 
source, gives an average value of (W(CO)) ~ 15 K kms -1 for 
an average velocity of Vo ~ _ 23 kms" 1 , which corresponds to 
a kinematic distance of either d — 2.1 or 7.7 kpc, assuming the 
kinematic velocity model of Wouterloot & Brand ( 1989), with 
the ambiguity arising from the galactic rotation curve. The 1° 
cloud diameter converts to densities of n ~ 56 cm 4 and n ~ 15 
cm -3 for the two distances, respectively, which translates to a 
total required cosmic ray energy of either ~ 3x 10 48 (for d = 2.1 
kpc), or, ~ 10 50 ergs for d = 7.7 kpc. The size of the y-ray 
source is determined by the size of the dense target (density 
profile) and the diffusion (propagation) length. The densities 
derived here will be required for the discussion below. 

(iii) The third CO solution involves the detection of the gi- 
ant molecular cloud (GMC) G303. 9-0.4 in the Carina arm at 
a distance of 12 kpc. The I - b CO map of G303.9-0.4 (in- 
tegrated over velocities in the 7 to 50 kms -1 interval) can be 
found in Grabelsky et al. (1988). Two adjacent CO maxima are 
visible, with the TeV source located on one of them, giving a 
TeV source averaged CO intensity of (W(CO)) ~ 8 K kms -1 . 
The target density at the contours corresponding to the TeV 
maximum is then n ~ 10 to 20 cm 4 , if a cloud scale size cor- 
responding to the FWHM of the integrated temperatures on the 
image is assumed. The required energy in cosmic rays is now 
Eqr ~ 3 x 10 50 ergs. Such energetics would require a relative 
powerful supernova explosion and it remains to be shown if an 
old SNR shell with 40 pc radius can still support cosmic ray 
production at the required level. 

An interesting counterpart for HESS J 1303-631 in connec- 
tion with CO solution (ii) may be the OB stellar association 
Cen OBI. Via the stellar winds of their member stars, OB 
associations are thought to supply considerable kinetic en- 
ergy into the surrounding media. Cen OB 1 contains ~20 stars 
of spectral type 09, suggesting that some member stars may 
have evolved beyond the supernova stage. The Wolf-Rayet 
star 0-Mus is also a member and may be a triple (WR + O 
+ O stars) system (Hartkopf et al. 1999 1. McClure-Griffiths 
et al. ( 120011 identify the Hi shell GSH 305+01-24 as being 
blown out by the OB association Cen OB 1 . The estimated dis- 
tance of GSH 305+01-24 at 2.2 kpc is consistent with that of 
Cen OBI (2.5 kpc, Humphreys 1978 1. Moreover, the estimated 
kinetic energy required to maintain the H i shell expansion at 
~ 10 x 10 51 ergs is met by that of Cen OBI at ~ 8 x 10 51 ergs, 
based on stellar wind luminosities and ages of member stars. 
Figure depicts the positional relationship between the H i 
shell, Cen OBI member stars and HESS J 1303-631. Assuming 
a canonical factor 0. 1 for kinetic energy to comsic-ray energy 
conversion, ~ 8 x 10 50 ergs would be available for cosmic -ray 
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production. This link between a TeV y-ray source and an OB 
association mirrors very closely that of the first unidentified 
TeV source TeV J2032+4130 and Cygnus OB2 discovered by 
the HEGRA collaboration (Aharonian et al. l2002all2"0"04at . 

Note, that the distances of Cen OBI and GSH 305+01-24 
coincide with the first kinematic CO solution of 2. 1 kpc (within 
uncertainties). A number of scenarios discussed below could 
give rise to the TeV emission: 

(a) The stellar winds from member stars of Cen OBI pro- 
vide conditions for the acceleration of cosmic rays, perhaps 
via shocks setup at the interfaces of the various stellar winds. 
These cosmic rays in turn could illuminate a nearby molec- 
ular cloud. If the average distance of HESS J1303-631 to the 
Cen OBI member stars is of order < 2°, the required energy 
to account for the TeV source, assuming all member stars con- 
tribute, would be 4 x 10 48 (2.0/0.32) 2 erg = 2 x 10 50 erg. This 
is consistent with the cosmic -ray energetics of Cen OBI dis- 
cussed earlier. Here, a source radius of 0.32° (95% or 2cr con- 
tainment radius) is assumed. A dense cloud situated at the lo- 
cation of HESS J 1303-631, could therefore be a valid solution. 
The energetics could also be met by a single or a few mem- 
ber stars, provided they were close enough. The CO data how- 
ever, also suggest the presence of numerous dense clouds in the 
vicinity of Cen OBI. The fact that TeV emission from just one 
location is seen could imply that HESS J 1303-631 is somehow 
special in relation to the distance to the local cosmic-ray accel- 
erator, and its geometry. Detailed modeling of the acceleration 
and diffusion of cosmic -rays in this region will be needed to 
verify this scenario. 

(b) As for situation (a) but the TeV emission may arise 
from a leptonic process such as inverse-Compton upscattering 
of local soft photons. The localised nature of HESS J 1303-631 
would then require some extra conditions namely: the soft pho- 
ton distribution (which may be dominated by UV photons from 
the member stars) is also localised; the acceleration and diffu- 
sion properties of electrons is such that the energetic electrons 
are confined. 

(c) If the oldest member stars already evolved into the su- 
pernova stage, the SNOB scenario of Montmerle ( 1979 1 may 
also explain the existence of HESS J 1303-631, but it remains 
then to be shown why this SNR is not seen at X-ray and radio 
wavelengths. 

5.4. Clumps of annihilating dark matter 

A possible source of TeV y-ray emission could also be a clump 
of annihilating dark matter particles (Moore et al. 11999 ) which 
would result in an energy spectrum with a cut-off at the mass of 
the dark matter particle. The galactic halo could contain a huge 
amount of such clumps comprised e.g. of neutralinos. 

Although not very likely, it is possible to reproduce in 
Monte Carlo simulations the observed flux with the annihila- 
tion radiation from high mass neutralinos, but it is not possible 
to reproduce the luminal profile measured from HESS J 1303- 
63 1 assuming a Navarro-Frenk-White (NFW) dark matter den- 
sity profile (Navarro et al.Q996 1. For the NFW and other com- 



monly used profiles a narrower luminal profile is expected than 
obtained from HESS J 1303-631. 



6. Summary & Conclusions 

A new unidentified TeV y-ray source HESS J1303-631 was 
serendipitously discovered in a dataset which was initially 
taken on the binary system PSRB1259-63/SS2883 which 
was also discovered at TeV energies (see a parallel paper: 
Aharonian et al. 2004c I. For the first time in TeV y-ray astron- 
omy, the detection and analysis of two sources within the same 
field of view is achieved which subsequently also occured in 
other pointings towards the galactic plane, showing the poten- 
tial of the new generation of ground-based experiments - such 
as H.E.S.S. - with the stereoscopic observation mode and its 
large field of view of ~ 5°. 

HESS J1303-631 was found to be clearly extended with 
a width of an assumed intrinsic Gaussian emission profile of 
cr = (0.16 + 0.02) °. The energy spectrum can be described by 
a power-law with a photon index of T = 2.44 + 0.05 stat ± 0.2 syst . 
The integral flux above 380 GeV was found to remain on a con- 
stant level of (17 ± 3) % of the flux from the Crab Nebula dur- 
ing the observations taken between February and June, 2004. 
Detailed studies on morphology and spatially-resolved energy 
spectra are underway. 

Up to now, no counterpart at other wavelengths was identi- 
fied. However, the location close to the galactic plane places 
HESS J1303-631 in the vicinity of a variety of possible ob- 
jects which might be involved in the production mechanisms 
explaining the observed TeV y-ray emission. A pulsar wind 
nebula powered by the nearby pulsar PSR J 130 1-6305 is one 
of the interesting candidates. Another possibility would be the 
Cen OBI stellar association at a distance of 2.5 kpc which 
provides enough kinetic energy by stellar winds. The interac- 
tion of an expanding supernova remnant shell with a molec- 
ular cloud as target material can also be seen as an interest- 
ing configuration for TeV y-ray production. Three molecular 
clouds are located along the line of sight of HESS J 1303-631: 
The Coalsack nebula at a distance of ~ 175 pc, a cloud in the 
Carina-Sagittarius arm at 2.1 kpc, resp. 7.7 kpc, and the giant 
molecular cloud GMC G303. 0-0.4 in the Carina Arm at a dis- 
tance of 12 kpc. Finally, the possibility of a clump of annihilat- 
ing dark matter particles was considered, which however can 
be excluded since the luminal profile of HESS J 1303-631 is not 
reproduced. 

In summary, no clear counterpart was identified to date. 
Therefore HESS J 1303-631 has to be considered as the second 
unidentified TeV source detected, following TEV J2032+4130 
discovered in the Cygnus region by HEGRA (Aharonian et al. 
2002al l2004ab . Meanwhile, the results of the H.E.S.S. galac- 
tic plane scan revealed more unidentified TeV y-ray sources 
(Aharonian et al. 2005b), thus further opening the door to 
a new class of (yet-unidentified) TeV y-ray sources. To fur- 
ther investigate possible production mechanisms and to under- 
stand this new region of the non-thermal universe future multi- 
wavelength observations (especially in X-rays) are essential 
and partially already initiated (Mukherjee & Halpern 2005). 
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